December 1968 


Brief 68- 10450 



AEC-NASA TECH BRIEF 



AEC-NASA Tech Briefs describe innovations resulting from the research and development program 
of the U.S. AEC or from AEC-NASA interagency efforts. They are issued to encourage commercial 
application. Tech Briefs are published by NASA and may be purchased, at 15 cents each, from 
the Clearinghouse for Federal Scientific and Technical Information, Springfield, Virginia 22151. 


Computer Program TRACK Performs Transient and/or Steady State 
Thermal Analysis with Coupled Fluid Flow and Heat Conduction 


The problem : 

To compute detailed transient and/or steady state 
fluid conditions (flow, fluid temperature, and pressure 
distributions) and spatial material temperature distri- 
butions for reactor components, and other types of heat 
exchange apparatus or components. The specified con- 
ditions are the geometric parameters of the flow system 
whiclvconsists of multiple, parallel fluid channels for 
cooling of the solid material, plena initial conditions, 
the solid body geometry which can be arbitrary in shape, 
and internal nuclear heat generation rates of the solid 
materials, if any. 


The solution: 

A computer program called TRACK was developed 
by combining a transient fluid flow computer code and 
the existing modified TOSS heat conduction code to 
perform the computation. 

How it's done: 

The transient solution at the end of a time step is 
obtained by iterating the channel wall temperatures 
between the fluid flow andthe heat conduction analysis. 
The procedure starts with a trial channel wall temper- 
ature distribution. The fluid flow calculation distributes 
the fluid between the various flow channels, calculates 
the convective heat transfer coefficients, the coolant 
temperature and pressure distributions along the chan- 
nels. Either the total flow rate or the system overall 
pressure drop between inlet and outlet plenum may be 
specified. The coolant temperatures and the heat trans- 
fer coefficients thus obtained are applied as boundary 
conditionsto the heat conduction computation. The re- 
sultant channel wall temperatures are compared to the 


trial channel wall temperatures used in the fluid flow 
calculation. Iterations continue until the channel wall 
temperatures are within a specified tolerance at the end 
of the time step. Computation then proceeds to the next 
time step. 

Similar procedures are used for steady state deter- 
mination. A transient problem can be started after a 
steady state solution is obtained. 

The program uses finite difference solutions to solve 
the governing transient fluid flow and heat conduction 
equations. It can handle up to 46 parallel flow channels 
with various sfngle phase fluids. The heat conduction 
model may be two or three dimensional and consist of 
different solid materials with temperature dependent 
properties. 

Notes: 

1. The program is written in Fortran IV for use on 
either the CDC-6600 computer or IBM-7094 with a 
64Kcore. 

2. Inquiries should be made to: 

COSMIC 
Computer Center 
University of Georgia 
Athens, Georgia 30601 
Reference: B68- 10450 

Patent status: 

No patent action is contemplated by AEC or NASA. 

Source: A. Y. Lee and M. D. Woods 
of Westinghouse Astronuclear Laboratory 
under contract to 
AEC-NASA Space Nuclear Propulsion Office 

(NUC-10189) 
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